Superconducting Instability in the Periodic Anderson Model 
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Employing a quantum Monte Carlo simulation we find a pairing instability in the normal state 
of the infinite dimensional periodic Anderson model. Superconductivity arises from a normal state 
in which the screening is protracted and which is clearly not a Fermi liquid. The phase diagram 
is reentrant reflecting competition between superconductivity and Fermi liquid formation. The 
estimated superconducting order parameter is even, but has nodes as a function of frequency. This 
opens the possibility of a temporal node and an effective order parameter composed of charge pairs 
and spin excitations. 



Introduction A number of Heavy Fermion materi- 
als display highly unusual superconductivity (Ott 1994, 
Grewe and Steglich 1991, Hess, Riseborough and Smith 
1993). It seems unlikely that conventional superconduc- 
tivity coexists with the strong local Coulomb correlations 
necessary to enhance the electronic mass in these sys- 
tems. Indeed, the specific heat jump at the transition 
scales with the normal state specific heat making it clear 
that pairing is between the heavy electrons. However, 
strong Coulomb correlations present no problem for un- 
conventional superconducting order parameters with ei- 
ther spatial (Sauls 1994) or temporal (Berezinski 1974, 
Balatsky and Abrahams 1992, Abrahams 1995) nodes. 
The latter can occur with either odd or even-frequency 
order parameters, which are interpreted in terms of com- 
posite condensates (Bonca and Balatsky 1994). Such an 
interpretation is supported by two sets of data: (i) power 
laws observed in physical properties below the supercon- 
ducting transition (Ott 1994, Sauls 1994), which contrast 
with the activated behavior of the conventional (node- 
less) s-wave order of, e.g., lead; (ii) the complex super- 
conducting phase diagrams of UPt3 and Ui-aTh^Beia 
(and possibly UBei3 itself, for which the penetration 
depth displays evidence of a secondary transition) (Ott 
1994, Sauls 1994). Finally, at least in UBei 3 (Ott 
1987) and CeCu2Si2 (Steglich 1996), the superconductiv- 
ity arises in a normal state which is clearly not a Fermi 
liquid. For T > T c in each of these materials, the lin- 
ear specific heat (C/T) rises with decreasing T while the 
resistivity decreases and has a high residual value at T c 
(typically 80-100 /if2-cm in the best samples of UBei.3). 

We have recently analyzed some static (Tahvildar- 
Zadeh, Jarrell and Freericks 1997a) and dynamic 
(Tahvildar-Zadeh, Jarrell and Freericks 1997b) proper- 
ties of two theoretical paradigms of heavy fermion materi- 
als: the single-impurity Anderson model (SIAM) and the 
periodic Anderson model (PAM). We found for metallic 
systems with an /-band filling nj«l and d-band filling 
n-d ^ 0.8, the Kondo scale for the PAM, To, is strongly 
suppressed compared to Tq IAM , the Kondo scale for a 
SIAM with the same model parameters. Consequently, 



the temperature dependence of the Kondo peak in the / 
spectral function is protracted in the PAM, i.e. it is much 
weaker than in the SIAM, consistent with experiments 
on single-crystalline Kondo lattice materials (Andrews 
1996). However the high temperature (T > T§ IAM ) 
properties of the two models are similar, so that Tq IAM 
is still the relevant scale for the onset of screening in 
the PAM whereas To, which bears no obvious relation 
to Tq IAM , is the scale for the onset of coherence where 
the moment screening is almost complete and the Fermi 
liquid begins to form. 

We interpreted our results as indicating the emergence 
of a single heavy band described by an effective Hub- 
bard model for the local screening clouds introduced by 
Nozieres (Nozicres 1985). Since the effective model is 
close to half filling and its hopping constant is strongly 
suppressed by the overlap of the screened and unscreened 
states of the f-electron moments, its Kondo scale, To, be- 
comes much less than Tg IAM . In this paper, we show 
that the regime of protracted screening is also associated 
with the occurrence of a pairing instability in the PAM. 

Formalism The PAM Hamiltonian on a £>-dimensional 
hypercubic lattice is, 

+ E (^4** + tffU*) + V J2 (4/<* + H - c -) 

+ E f/ ("/*t- 1 /2)(n/a-i/2) • (!) 

i 

where d(f)^J destroys (creates) a d(f) electron with spin 
a on site i. The hopping is restricted to the nearest 
neighbors and scaled as t*/2>/~D. U is the screened on- 
site Coulomb repulsion for the localized / states and V 
is the hybridization between d and / states. We choose 
t* as the unit of energy (t* w a few electron-volts, the 
typical band- width of conduction electrons in metals). 

We solve this model Hamiltonian with the dynami- 
cal mean field (DMF) method of Metzner and Vollhardt 
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(Metzner and Vollhardt 1989) who observed that the ir- 
reducible self-energy and vertex-functions become purely 
local as the coordination number of the lattice increases. 
As a consequence, the solution of an interacting lattice 
model in D = oo may be mapped onto the solution of 
a local correlated impurity coupled to a self-consistently 
determined host (Pruschke, Jarrell and Freericks 1995; 
Georges, Kotliar, Krauth and Rozenberg 1996). We 
employ the quantum Monte Carlo (QMC) algorithm of 
Hirsch and Fye (Hirsch and Fye 1989) to solve the re- 
maining impurity problem and calculate the imaginary 
time, local, one and two particle Green's functions. The 
maximum entropy method (Jarrell and Gubernatis 1996) 
can then be used to find the / and d density of states and 
the self-energy. 

We search for pairing between the f-electrons by cal- 
culating the f-electron pair-field susceptibility in the 
normal state. As shown in Fig. fjl this may be re- 
expressed in terms of the pairing matrix M follow- 
ing Owen and Scalapino (Owen and Scalapino 1971). 
M is formed from the irreducible superconducting 
(particle-particle, opposite-spin), local vertex function 
T sc , and the bare pair susceptibility x a (iu> n ,q) = 
Sj^ m ,k G '(^™' k ) G '(^™ ~ iw m ,q-k), with G(iw m ,k) 
the fully dressed lattice propagators; 



M = Vx°(iu>„,q)r£ m vW^Wq) ■ 



(2) 



The total pair susceptibility in terms of M is given by 
P SC (T) = VX°(*^ q) (1 - M)- 1 Vx°(ia, m , q) . (3) 

For the usual second order normal-superconducting tran- 
sition T c is obtained when the largest eigenvalue A of 
M reaches one, so that the pair susceptibility diverges. 
The corresponding eigenvector <fr(iu> n ) yields information 
about the superconducting order parameter. The pair- 
field susceptibility in this most singular channel can be 
projected out, using an appropriate frequency form fac- 
tor /„ = $(!w n )/ v /x n K,q), 



fr, 



(4) 



Results and Interpretation We introduce an effective 
hybridization strength T(u>), 



r(w) = Im £(w) 
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(•5) 



where G/(w) is the local / Green's function and 
is the local f-electron self-energy. T(oj) is a measure of 
the hybridization between the effective impurity in the 
DMF problem and its medium (for example, in the SIAM 
r(w) = irV 2 N,i(uj), where Nd(uj) is the d-band density of 



states). Fig. || shows the effective hybridization for the 
PAM near the Fermi surface /i = 0. In this figure the 
model parameters are chosen to be U = 1.5, V = 0.6, 
nt « 1 and three different d band fillings = 0.4 
(T = 0.014), n d = 0.6 (T = 0.054) and n d = 0.8 
(To = 0.16). A relatively small value of U/V 2 was chosen 
for presentation purposes; however, the features shown 
are also present for larger values of U/V 2 . (For the pa- 
rameter set of Fig. |^ we do not observe any supercon- 
ducting instability, although the hybridization functions 
look qualitatively very similar to the case where the in- 
stability exists). 

The feature we want to emphasize on is the dip in 
r(a;) at the Fermi energy which develops as the temper- 
atures is lowered. Since only the electronic states within 
about To of the Fermi surface participate in the screening 
(Nozieres 1985) this indicates a reduced number of avail- 
able screening states at the Fermi energy, and since the 
dip becomes narrower as rid — * 1, the effect is more dra- 
matic for small rid- Thus, as the temperature is lowered 
and rid < 0.9, the Kondo scale of the DMF effective im- 
purity problem is self-consistently suppressed and hence 
the "coherence" Kondo scale T is also suppressed. Con- 
comitant with this suppression, the Kondo length scale 
Ik ~ 1/To of the Kondo spin and charge correlations 
will increase dramatically. We have seen previously that 
these correlations can induce ferromagnetism in the PAM 
(Tahvildar-Zadeh et. al. 1997a) and superconductivity 
in the two-channel Kondo model (Jarrell, Pang and Cox 
1997). 

As discussed in the previous section, a second order 
superconducting transition is indicated when the largest 
eigenvalue of the pairing matrix exceeds unity. However, 
in our calculation, the largest positive eigenvalue does 
not approach unity for any value of U or temperature 
studied. Instead, for relatively large values of U > 2t* 
we see a different type of instability: At a high tem- 
perature (T > t*) we find generally that all eigenvalues 
are either small in magnitude or large and negative. As 
the temperature is lowered, the positive eigenvalues grow 
slowly; however, the most negative eigenvalue diverges 
at a temperature T*. Upon lowering the temperature 
further, the dominant eigenvalue (that with the largest 
absolute value) switches sign and becomes large (^ 1) 
and positive. It then decreases for some range of tem- 
peratures, but eventually increases again, diverges at a 
temperature T ; *, and switches back to large and nega- 
tive. The corresponding pair susceptibility P\(T) goes 
continuously through zero at the temperatures T*^ of 
the divergence of the dominant eigenvalue A. P\(T) is 
negative for T ; * < T < T* , indicating a pairing instability 
of the system in this channel. We find that these insta- 
bilities are degenerate throughout the Brillouin zone, i.e., 
P\(T — T*ns) vanishes for every q over the whole zone. 
We demonstrate this by plotting the corresponding local 
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pair susceptibility ph> cal (T) in Fig. || and observe that 
the instability prevails even at the local level. Hence this 
is a locally driven transition, consistent with a transi- 
tion driven by local dynamical correlations such as those 
responsible for Kondo screening. We suspect that the 
degeneracy in the Brillouin zone will be lifted in a finite 
D calculation by non-local dynamical correlations, such 
as spin waves, as has been suggested previously (Jarrell, 
Pang and Cox 1997). 

To understand how the vanishing of the pair suscep- 
tibility indicates a phase transition, remember that the 
inverse pair susceptibility is proportional to the curva- 
ture of the free energy /(A) as a function of the pair 
order parameter A, 1/P\(T) oc d 2 f '(A) / dA 2 . Thus, if 
Pa (T) < the normal state becomes thermodynamically 
unstable. The associated transition cannot be contin- 
uous, since this would require /(A) to become flat for 
small A (i.e. Pa (T) diverges) so that the order parameter 
may change continuously. Thus the observed transition 
is discontinuous. Furthermore, if P(T*) = 0, then T* is 
a lower bound to the transition temperature since when 
Px (T — T* = 0") = 0" the curvature of /(A) is diver- 
gent and negative, i.e. the free energy displays an up- 
ward cusp. This would compel the order parameter and 
the free energy to change discontinuously at T*, which 
involves an infinite energy at the transition. Hence, the 
actual transition occurs at a temperature T c > T* . Sim- 
ilar arguments can be made to show that T* is an upper 
bound to the discontinuous transition back to the normal 
metal. 

The dominant eigenvector of the pairing matrix at 
T = T* is even in frequency but has both positive and 
negative parts. To show the effect of the nodes on the 
order parameter in the ordered phase, we first introduce 
the Nambu-Gorkov matrix form for the local f-Green's 
function 



k lu n Z n l - e f T 3 - 0„ti - [iuJnI _ (ek+ed)T3] 

(6) 

where I, t\ and T3 are Pauli matrices. Z n = 1 — 
Im(Y,(iu> n )) /u) n is the wave function renormalization fac- 
tor and cj) n is the renormalized gap function. The f- 
electron order parameter is given by G/i2(t), the Fourier 
transform of the off-diagonal component of Gf(ico n ). In 
order to study this order parameter using our knowl- 
edge of the normal state, we assume that (f> n has the 
same temporal parity as the dominant eigenvector of 
the pairing matrix M right at T = T* (in the nor- 
mal phase). But we are not able to predict its actual 
value by this assumption. Hence, as the simplest ap- 
proximation, we assume that <p n = cf n where c is an 
unknown real normalization factor and /„ is the form 
factor constructed from the dominant eigenvector of the 
pairing matrix at T = T* . As c increases, G/i2(r) de- 



velops a suppression at r = (inset of Fig. ||). An or- 
der parameter with G/i2(0) « corresponds to a con- 
densate which excludes the simultaneous occupation of 
the same site by two electrons, minimizing its Hubbard 
energy. For the large values of U > 2t* where an in- 
stability is observed, we therefore assume that the con- 
densate with G/i2(0) « is the most stable one. For 
small r, an order parameter corresponding to an even- 
frequency <frn takes the form G/i2(t) » A + ^A 2 r 2 . If 
Ao = (/(t = + )/(r = 0)) = 0, then the superconduct- 
ing condensate is characterized by A 2 oc ([[f,H],H]f). 
In the strong-coupling limit, where the system can be 
mapped onto the Kondo lattice model, the commutators 
bring spin operators into the average (Bonca and Bal- 
atsky 1994), and if Ao = 0, the order parameter is effec- 
tively a composite of charge pairs and spin-excitations. 
Such an order parameter is compelling as it relates nat- 
urally to the complex phase diagram of the supercon- 
ducting heavy-fermion materials, in which the compe- 
tition and sometimes coexistence of antiferromagnetism 
and superconductivity is observed. However, it is unclear 
whether our estimate corresponds to the actual order pa- 
rameter of the broken symmetry state. A Monte Carlo 
simulation in the broken symmetry state is presently un- 
derway. 

Fig. ^ shows the resulting superconducting phase- 
diagram for the PAM with parameters U = 2, V = 0.5 
and nj « 1. The symbols correspond to the temperature 
where the pairing susceptibility P\(T) crosses zero and 
the largest diverging eigenvalue changes sign as discussed 
above. Note that the instabilities cease for > 0.9, 
leading to a phase diagram which suggests re-entrance 
into the normal state. Since T ; * is only an upper bound 
to T l c , an actual re-entrance transition need not occur, 
i.e. Tj could be vanishing. On the other hand, this 
would require the upper transition temperature (which is 
bounded from below by T*) to drop from a rather large 
value to zero upon a minute increase of n^. Although 
we cannot rule this out numerically, we believe this sce- 
nario to be unlikely. The upper part of Fig. ^ shows the 
Kondo scale for the PAM and SIAM with the same model 
parameters as stated above. We see that the supercon- 
ducting instability occurs only when To < Tq IAM , i.e., in 
the "protracted screening" regime. This again suggests 
that the pairing mechanism is related to underscreened 
moments and Kondo physics. Note that the upper tran- 
sition temperature is large, as is often the case for a mean 
field theory. We expect the non-local dynamical cor- 
relations of a finite dimensional system to significantly 
reduce our estimate for T*. Also, as noted previously 
(Tahvildar-Zadeh, Jarrell and Freericks, 1997) the pro- 
tracted screening regime and the associated supercon- 
ductivity vanish when the orbital f-degeneracy diverges. 
Thus, these phenomena are associated with finite orbital 
degeneracy. 

The shape of the lower phase boundary of the transi- 
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tion back to the normal metal is consistent with the rid 
dependence of the "coherent" Kondo scale To. This indi- 
cates that with the complete screening of the local mo- 
ments below T) the system either looses the mechanism 
driving superconductivity, or that the system can gain 
more free energy by forming a Fermi liquid. The shape of 
the upper phase boundary is more difficult to understand. 
T* increases upon lowering rid, in contrast to the simul- 
taneous decrease of the impurity Kondo scale Tq IAM . 
Hence, the onset of local moment screening alone (de- 
scribed roug hly by T^ IAM ) does not lead to supercon- 
ducting correlations. As observed by Nozieres (Nozieres 
1985), there are not enough states near the Fermi sur- 
face to completely screen all of the f-moments, especially 
as the d-band filling decreases. It is compelling to re- 
late T* to the temperature where the conduction band 
states available for screening are exhausted. Then, for 
T < T*, correlated lattice effects, such as superconduc- 
tivity or the magnetic polarons described by Nozieres, are 
required to quench the relatively large entropy associated 
with the unscreened moments. Better understanding of 
these correlations will be essential in a more quantitative 
explanation of the observed superconducting phase. 

Some understanding of these correlations is provided 
by examining the Matsubara frequency structure of the 
vertex function T sc (iu> n , iuj m ) which is minus the effective 
interaction (e.g. to lowest order in perturbation theory 
r sc (iu>„, iijJ m ) «—{/). Fig. U shows r sc corresponding to 
the point b in the phase diagram of Fig. ^ (b is just in- 
side the superconducting region). T sc (iuj n , iuj m ) is large 
and positive, indicating an attractive interaction, only 
at small frequencies, and it displays a central minimum. 
The width of the peak at each rid is roughly 2ttT* , which 
we have argued above to be related to Nozieres' "exhaus- 
tion" energy scale. Thus, the energy scale inferred from 
the susceptibility data is also present in the vertex func- 
tion. This again hints to Kondo/Nozieres screening as 
the mechanism driving the superconductivity. 

Fig. H shows the irreducible superconducting vertex 
function in the uj n — uj m direction (T sc (iu> n , at 
three different points of the phase diagram (cf. Fig. |J). 
We generally find that the vertex function is purely neg- 
ative for points outside the superconducting region such 
as c, whereas it has positive values at low frequencies for 
points inside the superconducting region such as a and b. 
This means that the effective interaction between parti- 
cles of opposite spin is attractive for the points inside the 
superconducting phase, but it is repulsive in the normal 
phase. Furthermore, the narrow width and the dip in the 
vertex function around zero frequency signifies that the 
attractive interaction is highly retarded. The central dip 
almost becomes a node at points well within the super- 
conducting phase (such as a) indicating that the static 
component of the attractive interaction is not significant. 
Therefore, the system would not gain energy by pairing 



two electrons on the same site at the same time. This is 
consistent with a node in the order parameter Gfriij) at 
t = as discussed above. 

Conclusions In conclusion, we have found a pairing 
instability in the infinite-dimensional periodic Anderson 
model which we interpret as a superconducting instabil- 
ity. This instability has several unusual features includ- 
ing: (i) A reentrant phase diagram reflecting competition 
between superconductivity and Fermi liquid formation; 
(ii) For most of the phase diagram, the upper transition 
temperature T* 3> To, indicating that superconductivity 
arises from a normal state which is clearly not a Fermi 
liquid; (iii) Although even in frequency, the dominant 
eigenvector of the pairing matrix has nodes, which opens 
the possibility of a composite order parameter with a 
temporal node; (iv) The superconducting instability oc- 
curs only in the protracted screening regime, where the 
effective hybridization diminishes and the Kondo length 
increases with decreasing temperature; (v) Within the 
superconducting regime, the vertex function displays a 
narrow feature with central node at low frequencies, con- 
sistent with a highly retarded interaction with a weak 
static component. These features, together with the fact 
that the instability occurs simultaneously over the whole 
zone (i.e., is locally driven), suggests that the pairing 
mechanism has its origins in the spin and charge corre- 
lations associated with Kondo/Nozieres screening of the 
f-electron spins. 
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FIG. 1. Owen and Scalapino's pairing matrix formalism 
(Owen and Scalapino 1971). The diagrams for the pair-field 
susceptibility (top) after integration over the internal mo- 
menta may be re-written (bottom) in terms of the pairing 
matrix M = \/^T sc y^ (Eq. 2). Here T ac is the irre- 
ducible two-particle self energy matrix and x° the bare pair- 
ing susceptibility. The square root is indicated by displaying 
only "half" of the electron Green's functions in the appropri- 
ate places. Second order phase transitions are signaled when 
the largest eigenvalue of M reaches unity. The corresponding 
eigenvector $ yields information about the superconducting 
order parameter. 




FIG. 2. Near-Fermi-energy (fi — 0) structure of hybridiza- 
tion function for the asymmetric PAM. The model parameters 
are U = 1.5, V = 0.6 and rif ~ 1. The dip at the Fermi-energy 
denotes a decrease in the number of states available for screen- 
ing, leading to a suppressed "coherence" Kondo scale To- The 
zero temperature hybridization for the single impurity model 
is shown for comparison. 
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FIG. 3. The local pair susceptibility of the dominant 
eigenvalue versus temperature for the PAM with (7 = 2, 
V = 0.5, n d = 0.6 and n f w 1. The susceptibility p^> cai (T) is 
positive in the temperature regimes T > T„ and T < T* . It 
goes smoothly through zero at T*^ and would be negative in 
between, indicating an instability in the corresponding pair- 
ing channel. The eigenvalue A diverges and changes sign at 
T*n (not shown). The corresponding eigenvector has both 
positive and negative parts. The inset shows the correspond- 
ing superconducting order parameter Gf\2(j) calculated at 
T* ~ 0.24 as described in the text. Since cf>„ has arbitrary 
normalization, we multiply it by a parameter c. As c in- 
creases, the order parameter develops a suppression at r = 0. 
This indicates that the two elements of the pair are not on the 
same site at the same time, minimizing the Hubbard energy 
of the condensate. 
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FIG. 4. Superconducting phase-diagram versus the d-band 
filling n d (bottom) for the PAM with U = 2, V = 0.5 and 
n/ ~ 1. The symbols denote the upper and lower bounds 
for the region where the leading eigenvalue of the pairing ma- 
trix M diverges (see text). The dashed line is a guide to the 
eye. The upper part shows the Kondo scales for the PAM and 
SIAM with the same parameters as for the lower part. Note 
that superconductivity only occurs where the lattice scale is 
suppressed relative to the impurity scale, and that the up- 
per instability temperature T£ > T$ AM , indicating that the 
superconductivity emerges from a state which is not a Fermi 
liquid. 




FIG. 5. The superconducting vertex function correspond- 
ing to the point b of Fig. ^. The z axis is in units of the 
function value at the central minimum. Note that the ver- 
tex function is large only at low frequencies, with a width 
u n f» 2-kT*. 
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FIG. 6. The superconducting vertex function in cu n — aim 
direction at three different points of the phase diagram (cf. 
Fig. ^]). The vertex function is purely negative for point c out- 
side the superconducting region whereas it has positive values 
at low frequencies for points a and b inside the superconduct- 
ing region, illustrating the attractive pair interaction in this 
region. Note that the narrow peak width and the dip around 
zero frequency implies strong retardation of the interaction. 



